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ABSTRACT 

A well-tested and validated Gibbs sampling code, that performs component separation and CMB 
power spectrum estimation, was applied to the WMAP 5-yr data. Using a simple model consisting 
of CMB, noise, monopoles and dipoles, a "per pixel" low-frequency power-law (fitting for both am- 
plitude and spectral index), and a thermal dust template with fixed spectral index, we found that 
the low-€ (€ < 50) CMB power spectrum is in good agreement with the published WMAP5 results. 
Residual monopoles and dipoles were found to be small (< 3 /iK) or negligible in the 5-yr data. We 
comprehensively tested the assumptions that were made about the foregrounds (e.g. dust spectral 
index, power-law spectral index prior, templates), and found that the CMB power spectrum was in- 
sensitive to these choices. We confirm the asymmetry of power between the north and south ecliptic 
hemispheres, which appears to be robust against foreground modeling. The map of low frequency 
spectral indices indicates a steeper spectrum on average (/3 = —2.97 ±0.21) relative to those found at 
low (~GIIz) frequencies. 

Subject headings: cosmic microwave background — cosmology: observations — radio continuum: ISM 



1. INTRODUCTION 

Accurate separation of diffuse Galactic foreground 
emissions is one of the major challenges of cosmic 
microwave background (CMB) data analysis. Fore- 
ground emission must be subtracted accurately and er- 
rors should be correctly propagated through the data 
analysis pipeline to extract precise cosmological infor- 
mation. Furthermore, the use of Galactic masks re- 
duces the sky area that can be exploited for cosmo- 
logical analyses thus increasing the sample variance at 
large angular scales. At angular scales > 1°, the domi- 
nant foregrounds are from synchrotron, free-free and vi- 
brational dust emission from the Galaxy. In addition, 
there is evidence for an additional "anomalous" fore- 
ground at low freque ncies (~ 10 — 60 GH z), which is still 
not well understood dDavies et al. 2006; Dick inson et al. 
2006, 2007, 2009: Bonaldi et al. 2007; Hildcbrandt et al. 
20071 ,Dobler fc Finkbeiner ,2008a: .Gold et al...2009l l 
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B ased on earlier p io neerin g work bvlJewell et al] (l2004f) 
and 'Wan delt et all ('2004'), 'Eri ksen et all (|2008af ) de- 
scribe a Gibbs sampling algorithm that can perform com- 
ponent separation and power spectrum estimate at the 
same time, by sampling from the joint posterior proba- 
bility distribution for all fitted parameters i.e. it is based 
on Bayesian statistics. This method has the advantage of 
full propagation of error information through to the CMB 
power spectrum (errors for each estimate) but also al- 
lows considerable flexibility in fitting of foreground model 
parameters. These can include fixed index templates, 
free templates at each frequency, monopoles/dipoles, and 
any spectral function (e.g. power-law at each pixel). 
A detailed description of the algorithm and the imple- 
mentation of the parall elized Fortran90 code ( henceforth. 
Commander) is given bv lEriksen et all (|2008a[ ). The code 
also takes into accounts finite bandpass responses at each 
frequency. 

In a previous work (jEriksen et al.l l2008b). 
Comma nder was use d ana lyze the 3-year WMAP 
data ijHinshaw et al.l l2007f ). at a resolution of 3° 
(Nside = 64). The analysis confirmed the accuracy of 
the CMB power spectrum up to = 30 and provided 
estimates of the foreground components including a map 
of the low frequency foreground amplitude and spectral 
index at each pixel. Due to the limited number of 
frequency channels and frequency range of the WMAP 
data, a number of assumptions had to be made such 
as using a thermal dust template with fixed spectral 
index and a prior on the foreground spectral index. 
Furthermore, they discovered the existence of significant 
monopoles (offsets) and dipoles in the 3-year WMAP 
data. Although this did not significantly affect the 
cosmological interpretation (but did affect the details of 
the foregrounds) it was later confi rmed to be an error i n 
the processing of the 3-year data (jHinshaw et al.ll2009( l. 

In this paper we analyze the new WMAP 5-year data 
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using the Commander code and confirm the CMB resuhs 
from the WMAP team. ij2] outUnes the data and pre- 
processing along with a brief summary of diffuse fore- 
ground components and the templates and models that 
are used to fit them. The basic WMAP5 analysis is pre- 
sented in f|3l In §4] we investigate the sensitivity of the 
results on a variety of different foreground models. We 
show that the results are largely insensitive to the as- 
sumptions made. A comparison of \ow-£ CMB and fore- 
ground spectral index maps is made in fJ51 Conclusions 
are given in ^ 

2. DATA 

2.1. WMAP data and preprocessing steps 

The main input data are the 5-yr WMAP tempera- 
ture maps provided on Lambd a^^ . We perform the s ame 
pre-processing steps taken bv lEriksen et all ()2008bf ). to 
allow the Gibbs sampling to be performed relatively effi- 
ciently. The 1° smoothed maps frequency maps were fur- 
ther smoothed to a common resolution of 3° and down- 
graded to a HEALPix^^ resolution of Ngidc = 64, corre- 
sponding to a pixel size of 55'. The single frequency band 
maps are at nominal frequencies of 23, 33, 41, 61 and 
94 GHz. We assume a top-hat bandpass response func- 
ti on with width s given by the radiometer specifications 
in lJarosik et al.l (|2003f ). An additional 2 /iK white noise 
is added to the maps to regularize the noise covariance 
matrix at high ^-values and thus reduce the number of 
conjugate gradient descent (CGD) iterations when solv- 
ing the map-making equation. The additional noise has 
a negligible effect for £ ^ 50. 

Note that the amount of smoothing and additional 
noise is somewhat arbitrary, and was a choice that 
was made to balance the accuracy of the results and 
the amount of CPU time. A typical analysis takes 
~ 1 minute per Gibbs iteration on a modern 3 GHz CPU. 
Thus for a single run, with 10000 Gibbs samples dis- 
tributed over 4 chains (4 x 5 = 20 C PUs in total' ), takes 
a total of ~ 40 hours. We note that I Jewell et all ()2008D 
have recently developed a Markov Chain Monte Carlo 
(MCMC) algorithm that effectively solves the problems 
associated with direct Gibbs sampling in the low signal- 
to-noise regime. 

2.2. Foreground components and templates 

Diffuse foreground emission consists of at least 3 well- 
known components: synchrotron, free-free and thermal 
(vibrational) dust. In additional to this, there is signifi- 
cant evidence for an anomalous component, possibly due 
to spinning dust. We briefly review the components and 
simple models that are used when fitting to real data. 

Synchrotron radiation is emitted by relativistic cosmic 
ray electrons spiralling in the Galactic magnetic field. 
If the energy spectrum of the cosmic rays is a power- 
law, dN/dE = NqE^p, then the intensity is given by 

= LNoB^ffj'^^^^^i^-(P-^^/^a{p), where a{iy) is a weak 
function of frequency, L is the length along the line-of- 
sigh of the emitting volume, BcS is the effective mag- 
netic field strength and v is the frequency. A power-law, 
T(i/) = Aiy^", is therefore a reasonable approximation 

http://lambda.gsfc.nasa.gov/ 
http://healpix.jpl.nasa.gov 



for synchrotron radiation with /3s ^ —2.7 at radio wave- 
lengths. At higher frequencies, a single source of elec- 
trons will lose energy faster (spectral aging) resulting in 
a steepening of the spectral index, by ~ 0.5. However, 
multiple synchrotron components along the same line-of- 
sight will tend to result in a flattening of the spectral 
index at higher frequencies. At WMAP frequencies, the 
values of Ps ~ —3.0 ± 0.3 have been reported. The steep 
spectral index makes low frequency radio maps ideal tem- 
plates for synchrotron radiation. Th e most widely used 
is the all-sky 408 MHz map of Hasl am et al.l ()1982D . 

Free-free (bremsstrahlung) emission is emitted by 
free electrons in warm (Tg ss 10'' K) gas. At ra- 
dio frequencies, the radio spectrum has the form Tff oc 
j.-i.35j^-2.i^^ where EM is the Emission Measure. 
The spectrum therefore closely follows a power-law with 
a spectral index of /3ff = — 2.1. A more accurate cal- 
culation, taking into account the Gaunt factor, gives 
a slightly steeper index at WMAP frequencies, with 
Pff ~ —2.15 and some weak depende nce on Te an d other 
environmental conditions (Dickinso n et al.ll2003[ ). How- 
ever, without a large lever arm in frequency, the separa- 
tion of synchrotron and free- free emission can be difficult. 
An alternative method is to use a spatial template of the 
free-free. Recombination lines (e.g. Ha) are also pro- 
portional to EM and therefore provide a reliable tracer 
of free-free emission, if Te is known; at high latitudes, 
Te « 8000 K. The major limitation when using optical 
lines is the visual dust extinction. Fortunately, at high 
Galactic latitudes, the extinction by dust is small. Full- 
sky maps a re now available including correctio ns for dust 
extinction (|FinkbeineHl2003t [Dickinson et~an i2003V 

Thermal (vibrational) dust emission arises from the 
heating of dust grains by the UV interstellar radiation 
field. At an equilibrium temperature of Td « 18 K, it 
emits mostly in the far infra-red, peaking at A « 100 /im. 
For CMB experiments, it is only a major foreground at 
v > 100 GHz. A common approximation to the ther- 
mal dust spectrum is a modified black-body: T{iy) oc 
v'^B{v, Td) where e is the emissivity index, with typical 
values - 2 (Table [I]). At WM^P frequencies, we see the 
Ray leigh- Jeans tail of the thermal dust, and hence the 
spectrum can be approximated by a power-law, based on 
the emissivity index, e. However, thermal dust is rela- 
tively weak at WMAP frequencies, therefore spatial tem- 
plates are often used as a tracer instead. The most widely 
used t emplate is the dust prediction of Finkb einer et al.l 
(jl999( l. relative to 94 GHz. Based on IRAS/DIRBE data, 
it represents a global fit to the all-sky data, with 2 dust 
components with different temperatures and emissivities. 

Anomalous emission is still not well understood. It ap- 
pears to be an excess of emission at ~ 10 — 60 GHz and is 
tightly correlated to FIR dust emission. The best expla- 
nation is in terms of electro-dipole radiation from rapidly 
spinning small dust grains. Theoretical models predict 
a peaked (convex) spectru m with a peak frequency in 
the range ^ 20 - 4 GHz (iDraine fc Lazarian|[T998allbl : 
lAli-Haimoud et al.l2008 1 . The anomalous emission is dif- 
ficult to separate spectrally in the WMAP channels since 
it is falling off quickly in the lower WMAP bands, with a 
spectral similar to synchrotron. As with free-free, a sep- 
aration based on spatial morphology has been the most 
successful way to separate the anomalous emission with 
the other components, due to the remarkably tight cor- 
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relation with FIR templates (jPavies et al.ll200^ . 

In this paper, we use standard foreground templates to 
trace diffuse foreground components (also downloaded 
from the Lambda website). For the synchrotron emis- 
sion, we use the dcstripcd and sour ce-removed removed 
NCSA version of the Haslam et all ([l982) 408 MHz all- 
sky map. For free- free emission, we use the all-sk y 
dust-corrected Ha map produced by iFinkbeineil (|2003D . 
We check the consiste ncy with the similar template by 
iDickinson et al.l (|2003) . For thermal dust, w e use the ex- 
trapolation of model 8 of Finkbeiner et all ((1999.) (here- 
after FDS99) to 94 GHz. We check the consistency with 
the original 100 /im intensity map, and also the dust 
reddening E{B — V) map. All maps were smoothed to a 
common 3° resolution and downgraded to iVsidc = 64. 

3. BASIC WMAP5 ANALYSIS 

3.1. Signal model and priors 

To begin w i th, we apply the signal model as used by 
lEriksen et all (|2008bD to the WMAP 5-yr data. For this 
case, the signal model being fitted to the data, is of the 
form:- 



T^{p) = s{p) +ml+J2 ■ n(p)] 



1=1 



+b [t{p)a{,^) {iy/4''^'r\ + f{p)a{iy) {vlv^f (p) 

The thermodynamic temperature Tu{p), in each pixel 
p, consists of terms from the CMB (assumed to be Gaus- 
sian and defined by a power spectrum up to £ = 150), 
monopoles and dipoles at each frequency, a single dust 
template amplitude, 6, based on the FDS99 model 8 
prediction t(p) and with a fixed index of /3d = +1-7, 
and a single power-law (amplitude and spectral index) 
for each pixel in the R-J convention; a{v) is the con- 
version from thermodynamic to R-J temperature units. 
The power-law accounts for all the low frequency fore- 
ground emission. We adopt a weak uniform prior on the 
power-law spectral index of —4.0 < /3 < — 1.0 and a more 
constraining Gaussian prior with mean p = —3.0 and 
standard deviation A/3 = 0.3. These are multiplied by 
the appropriate Jeffreys prior to account for the different 
volume of likelihood space for different spectral indices 
(|Eriksen et al.|[2008aD . The mask is a downgraded ver- 
sion of the standard Kp2 mask (not including the major- 
ity of extragalactic sources) where any masked sub-pixels 
result in a masked pixel at iVgido = 64, resulting in 14.4% 
of pixels being masked. 

3.2. CMB results 

The CMB temperature power spectrum from 
Commander and the WMAP team is shown in Fig. [TJ The 
maximum likelihood value and 68% confidence regions 
were calcula ted using the " Gaussianized Blackwell-Rao" 
estimator of iRudiord et al.1 (2008). The power spectrum 
is in very good agreement with those from WMAP 3-yr 
data and those from the WMAP team fjNolta et al.l 
[2009I) . on a i-hy-l basis. The error bars are not directly 
comparable since we calculate the 68% confidence limit 
for each Cf while we plot just the Gaussian symmetric 
error bars from the WMAP team. 

To determine the impact on the cosmological parame- 
ters we ran the derived WMAP5 power spectrum through 



COSMOMC^"' ([Lewis &: Bridld I2002D for the s tandard 6- 
parameter ACDM model ( Dunklev et al.ll2009D. In asim- 
ilar way to the WMAP team, and also lChu et all (|2005D . 
we use the Blackwell-Rao estimates for i < 32, and fix 
the Cf's for £ = 33 — 50 during the run, corresponding 
to a constant in the log-likelihood. The Commander and 
WMAP likelihoods for the six-parameter set are com- 
pared in Fig. [2l The distributions are almost identi- 
cal. The Commander marginal means and standard devia- 
tions for these parameters are fib/j^ = 0.02285 ±0.00063, 
r2cdm/i^ = 0.1089 ± 0.0064, as = 0.793 ± 0.037, Ha = 
72.5 ± 2.7, = 0.9675 ±0.0146 and r = 0.0898 ±0.0177. 
These are in excellent agreement with the WMAP 5-yr 
Hkelihood values: fib/i^ = 0.02278 ± 0.00062, l^cdm^^ '= 
0.1093 ± 0.0067, erg = 0.793 ± 0.038, Hq = 72.2 ± 2.7, 
n, = 0.9645 ± 0.0147 and t = 0.0891 ± 0.0175. 



3.3. Foreground results 

The FDS99 dust template amplitude is 0.949 ± 0.003 
relative to 94 GHz. This is in good agreement with the 
WMAP 3-yr results and indicates that the dust model 
of FDS99 is a remarkably accurate prediction when av- 
eraging over the high latitude sky. Note that the formal 
error is likely to be under-estimated, since it is based 
on the contribution from the noise model (2 /iK white 
noise) and does not include model-dependent uncertain- 
ties. The CMB map and low frequency foreground am- 
plitude and power-law maps are shown in Fig. [3] By 
eye, the maps look similar to those produced in the 3-yr 
analysis. The error in the CMB map is dominated by the 
dipole aligned with the Galactic mask (dipole X) while 
the foreground errors increase at lower latitudes due to 
the increased foreground levels. 

A recent improvement in the Commander code is the 
ability to fit for foreground parameters inside the Galac- 
tic mask, without affecting the CMB power spectrum. In 
masked regions, the large-scale CMB power and any tem- 
plates (foregrounds, monopoles/dipoles) are subtracted 
first, based on the fits of the unmasked regions. The 
residual is then fitted by any remaining foreground pa- 
rameters on a pixel-by-pixel basis. In this case, the low 
frequency power-law (amplitude and spectral index) is 
fitted to each pixel inside the mask. The result is a full- 
sky map of the foregrounds (Fig. Except for a small 
discontinuity between at the mask boundary, the spec- 
tral index map is within the expected range. The dom- 
inance of free-free emission, with a theoretical spectral 
index /3g « —2.1, is apparent in the Galactic plane and 
in well-known star- forming regions such as Orion ((/, 6) ^ 
(210°, -10°)), the Gum nebula {{I, b) ~ (260°, -10°) and 
Ophiucus {{l,b) ~ (5°, +25°). The average spectral in- 
dex for the whole sky is /3 = -2.90 ± 0.26. At high 
latitudes (|6| > 20°), the spectral index has a range of 
-3.8 < f3 < -2.1 with an average of /? = -2.97 ± 0.21. 
The average is therefore very close to the prior mean but 
with a standard deviation significantly less, suggesting 
that the actual variation in /? is less than the standard 
deviation of the prior, A/3 = 0.3. In H4.4.2I we will ex- 
plicitly test the sensitivity of the results on the prior. 

3.4. Monopoles and Dipoles 
http:/ /cosmologist.info/cosmomc/ 
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Signal components and models. 



TABLE 1 

Spectra are given in terms of a Rayleigh-Jeans temperature, Trj. 



Component 



Models 



CMB 



1. Gaussian random field with known frequency spectrum. 

2. Uncorrelated pixels, CMB spectrum: Trj = Ax (e'"'/'=T 



Synchrotron 



1. Spatial template: Haslam et al. 408 MHz all-sky map 

2. Spatial template: K — Ka WMAP difference map 

3. Uncorrelated pixels, power-law: Tpj^j = Ax (iz/vq)^' 

4. Uncorrelated pixels, power-law with running: Tj^j = Ax (!y/;/Q)A+C'log(''/i'o) 
Free- free 1. Spatial template: Ha maps (Finkbeiner ^003; ^i^k^ 

2. Uncorrelated pixels, fixed index power-law: T^j = A X {v/vq)~'^-^^ 
Anomalous dust 1. Spatial template: Dust map extrapolated to 94 GHz ^Finkbeiner_et_aL ^999) 

2. Spatial template: K — Ka WMAP difference map 

3. Uncorrelated pixels, power-law: Trj = Ax {v/vq)I^ 

Thermal dust 1. Spatial template: Dust map extrapolated to 94 GHz (Finkbeiner_et_al^^99^ 

2. Uncorrelated pixels, power-law: Trj = Ax {v/vo)^'^ 

3. Uncorrelated pixels, modified black-body Trj = Ax i/'' B{i/,T^) 
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Fig. 1. — Top: CMB temperature power spectrum obtained by Comma nder (red filled squ ares) and the WM^P team (black filled circles) 
based on the WMAP 5-yr data. The best-fit ACDM model spectrum of INolta et al.l l| 200i|) is shown as a dotted line. The error bars show 
asymmetric 68% confidence limits for Commander, and Gaussian Icr uncertainties for the WMAP spectrum. Bottom: Difference between the 
WMAP team and Commander power spectra. 



An important resu lt of the WMAP 3-yr analysis of 
lEriksen et all (l2008br ) was the discovery of significant 
residual monopoles and dipoles in WMAP maps. This 
was later confirmed to be due to an error in the pro- 
cessing of the maps when estimating the offset in each 
map base d on a model of the foregrounds at the Galac- 
tic poles (jHinshaw et al.ll2009f ). The best-fit monopoles 
and dipoles in the WMAP 5-yr data are given in Ta- 
ble [51 The magnitude of the monopoles and dipoles is 



clearly smaller than those found in WMAP 3-yr data 
indicating the improved estimation of these terms. How- 
ever, there still may be small monopole/dipole residuals 
in the WMAP 5-yr data. In particular, the monopoles 
appear to be positively biased at ~ 3 /xK in four of the 
five bands, while at Q-band the monopole is 5.7±0.5 /iK. 
The errors are purely statistical, based on the 2fiK white 
noise added to the maps, but there will be additional 
sources of error from unaccounted instrumental noise 
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Fig. 2. — Comparison of marginal probability distributions from WMAP 5-yr data of Commander (solid line) and the WMAP tca,m (dashed 
line) for a six-parameter cosmological model. 




Fig. 3. — Marginal mean maps (top) and r.m.s. (bottom) maps in Galactic coordinates. From left to right, the panels show the CMB, 
the foreground amplitude and the foreground spectral index, respectively. The signal model and Galactic mask is the same as used by 
lErikscn ct al. (2008b); see text for details. The Kp2-based mask is indicated in the grey shaded region of the CMB map. The color scale 
in the CMB maps is linear in units of /iK while the foreground amplitude is histogram equalized and in units of /iK normalized to 23 GHz. 
The new version of the Commander code allows the foreground parameters to be sampled inside the mask. 



and more importantly, modeling errors. It is possible 
for example, for power to alias from the foregrounds to 
the monopoles/dipoles and vice versa. We reduce this 
effect by applying an additional constraint on the or- 
th ogonality between the se two components as described 
in lEriksen et al.l ()2008af ). However, with such a simple 
model of a single power-law at each pixel, which accounts 
for several Galactic components (synchrotron, free-free 
and possibly anomalous dust emission) , there will be ad- 
ditional modeling errors. We will test the stability of 
these quantities by making adjustments to the model 
(masks, priors, templates etc.) and observing the spread 
in these values. 

To see how important the monopole and dipole terms 
are, we repeated the above analysis but not fitting for 
any monopoles and dipoles i.e. assume they are zero. 
The most obvious effect is in the foreground spectral in- 
dex map, shown in Fig. [4l Besides for differences inside 



TABLE 2 

Mean and r.m.s. of monopole and dipole posteriors in the 
WMAP 5-YR MAPS. Top: For the basic model using 1 "index" 
TEMPLATE (sEE Q. Bottom: For the sa me model but w ith 2 
"index" templates, including t he FinkbeinerI II2003I1 Hq 

TEMPLATE (SEE i|4.2.1ll . 



Band 


Monopole 


Dipole X 


Dipole Y 


Dipole Z 




(/.K) 


(mK) 


(/.K) 




K-band 


3.2 ±0.5 


2.3 ± 1.2 


0.4 ± 0.8 


-1.1±0.1 


Ka-band 


3.6 ±0.5 


1.0 ± 1.2 


-2.6 ±0.8 


0.5 ±0.1 


Q-band 


5.7 ±0.5 


1.9 ± 1.2 


-2.1 ± 0.8 


0.8 ±0.1 


V-band 


2.6 ±0.5 


2.4 ± 1.2 


0.6 ±0.8 


-1.3 ±0.1 


W-band 


2.6 ±0.5 


2.1 ± 1.2 


0.9 ±0.8 


-1.6 ±0.1 


K-band 


0.5 ±0.5 


2.9 ± 1.2 


1.2 ±0.8 


-0.8 ±0.1 


Ka-band 


0.3 ±0.5 


2.0 ± 1.2 


-1.3 ±0.8 


0.8 ±0.1 


Q-band 


2.6 ±0.5 


2.8 ± 1.2 


-0.7 ± 0.8 


1.1±0.1 


V-band 


-0.0 ±0.5 


3.0 ± 1.2 


1.6 ±0.8 


-1.0 ±0.1 


W-band 


-0.0 ±0.5 


2.7 ± 1.2 


1.6 ±0.8 


-1.4 ±0.1 
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the masked region (particularly in the Galactic center re- 
gion), there is a large-scale gradient running across the 
map, which is unlikely to be real, and is just a conse- 
quence of residual monopoles and dipoles in the WMAP 
5-yr maps. To verify this conclusion a WMAP simulation 
was made, with and without the best-fitting monopoles 
and dipoles. The spectral index map was seen to have 
a similar gradient. This is also obser ved to some degree 
in the spectral index maps derived bv lGold et al.l ()2009D 
using a MCMC algorithm, presumably a result of not 
accounting for residual monopoles / dipoles (see ii5.3l and 
Fig. [HI). 



Chi-squarcd 




Foreground index mean 




Fig. 4. — Low- frequency foreground spectral index mean map 
for WMAP 5-yr data when not fitting for monopoles and dipoles. 



4. WMAPb ANALYSIS WITH VARYING MODELS AND 
ASSUMPTIONS 

4.1. Galactic masks 

The Galactic mask is extremely important for obtain- 
ing the optimal CMB estimates by masking out bright re- 
gions dominated by foreground emission (which becomes 
increasingly complex, particularly at low Galactic lati- 
tudes) but still allowing the maximum sky coverage (to 
reduce cosmic variance) . From the analysis with the Kp2 
mask, the map gives a handle on how good the signal 
model and where it is a poor fit. Fig. [5] shows the map of 
X^. Most pixels have < 15, which indicates a model 
consistent with the data and the errors; statistically, a 
value of 15 corresponds to a model that is rejected at 
the 99% confidence level. There is a small increase in 
along the ecliptic plane, due to the additional unmodeled 
instrumental noise, although this is an improvement over 
the WMAP 3-yr results due to the additional data. A 
few pixels have particularly large values (x^ > 50). 
An obvious mask choice would be to only include pixels 
with a x^ below some limit. We tried a mask for x^ < 30 
and found that the CMB power spectrum remained un- 
changed. 

We explore the sensitivity of the derived CMB and 
foreground parameters by first using the exact same sig- 
nal model as in Sj3] but employing five different masks. 
The masks are i) Kp2 mask (14.4%) as used in ^ ii) 
a 10° smoothed Kp2 mask (31.4%), iii) WMAP5 KQ85 
temperature mask (17.5%), iv) WMAP5 temperature 
processing mask (6.5%) and v) no mask (full-sky). The 
coverage of each mask is shown in Fig. [6] 



Fig. 5. — Map of for the basic WMAP 5-yr analysis. Most 
pixels are at < 15 and are consistent with the model, given the 
assumed errors. There is a slight increase in the x^ in the ecliptic 
plane due to the additional unmodeled instrumental noise. A few 
pixels close to the Galactic plane have > 50i with one pixel at 



138. 





Fig. 6. — Map showing the coverage of Galactic masks used. The 
masks are Kp2 (light blue), a 10° smoothed Kp2 (orange), WMAP5 
KQ85 temperature mask (green) and the WMAP5 temperature 
processing mask (dark blue). 



Fig. [7] shows the CMB temperature power spectra for 
each mask. The signal model is the same in each case. 
For full-sky (no mask) there are significant foreground 
contributions across the entire €-range considered here 
{£ — 2 — 50) with an excess power of w 500 /iK^. With the 
WMAP processing mask (6.5% masked out), the spec- 
trum is quite close to the spectra employing larger masks. 
However, there is some positive bias at some ^-values 
(e.g. ^ = 22 — 27), indicating that for this model, this 
mask is not adequate to mask out the Galaxy and thus 
have a negligible effect on the power spectrum. The tem- 
perature KQ85 mask and Kp2 mask power spectra are 
in excellent agreement. It is interesting to see that the 
largest mask, the smoothed Kp2 mask, produces a power 
spectrum that is slightly lower than the others, particu- 
larly at a few specific values (e.g. ^ = 10, £ = 25 — 28). 
This is partly attributed to the approximate G aussian- 
ized Blackwell-Rao estimator of iRudiord et al.l pOOS) . 
due to the extra correlations introduced by larger masks 
(> 30% cut). Even so, inspection of the marginal Ci his- 
tograms does confirm that a few specific ^-values (par- 
ticularly £ = 10) are lower, but no more than 2(t, than 
for the larger mask, and are likely to be due to residual 
foreground contamination at low Galactic latitudes. We 
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checked that the impact on the cosniological parameters 
is neghgible (<< la for all parameters). 

4.2. Ha templates 

A major issue with the signal model used so far is 
that the majority of the foreground emission is fitted 
by a single power-law at each pixel. This is a broad as- 
sumption given that we know that there are at least two 
distinct Galactic components emitting at the lower fre- 
quencies: synchrotron radiation (which has a spatially 
varying spectral index) and free-free emission (which as 
a flatter spectral index but varies little with position). 
There is also substantial evidence for an anomalous dust- 
correlated components, possibly due to spinning dust 
grains, which is strong in the lower WMAP bands. The 
anomalous component appears to ha ve a spectral index 
of /3 « -3.1 between 20 and 40 GHz (jPavies et al.ll2006f) 
thus one might expect a single power-law to be ade- 
quate for both components. However, the free-free emis- 
sion has a spectral index « —2.15 at these frequen- 
cies ([bickin son et al. _2003) . If the free- free emission is 
significant in the higher WMAP bands, the power-law 
model will clearly not be a good model and the CMB 
estimates could be affected. Fortunately, there are fuU- 
sky free-free te mplates available, based on Ha surveys 
(jPickinson et a l. 2003; Finkbeincr 2003). These are ex- 
pected to be good tracers of the free-free emission for 
most of the sky, with the exception of low latitudes where 
dust absorption will affect the Ha intensity. 

4.2.1. Ha fixed index template 

We included a free-free teniplate based on Ha 
(|Finkbeineij|2003HDickinson et al.ll2003l ) as an additional 
template with a fixed spectral in dex i3s = ~2-15. This 
was not included in the analysis of'Erikse n et al.l (|2008br ) 
because simulations showed that with this model, for the 
WMAP frequency range, the code did not reproduce the 
exact input values. In particular, the free-free amplitude 
was underestimated by « 10% (depending on the typi- 
cal assumed synchrotron spectral indices) due to aliasing 
between the power-law component and the free-free tem- 
plate. This is due to a degeneracy between synchrotron 
and free-free emission over the frequencies: With 

effectively only three available frequencies (K,Ka,Q) and 
two free components (synchrotron and free- free), the re- 
sulting joint posterior becomes strongly non-Gaussian 
and asymmetric. As a result, the corresponding univari- 
ate marginal means are no longer unbiased estimators, 
and not useful as point e stimators. Fo r an ex plicit ex- 
ample of this effect, see lEriksen et al.l (|2007[ ). In the 
simulation, the other components (dust template ampli- 
tude, monopoles and dipoles) were reproduced correctly 
and more importantly, the CMB power spectra estimates 
were still consistent within a small fraction of the formal 
errors, giving us confidence that this could still improve 
the foreground subtraction. Additional bands, particu- 
larly at lower frequencies, can improve this situation. 

When applied to WMAP 5-yr data, the foreground 
template amphtudes were 0.998 ± 0.003 for the FDS99 
dust template at 94 GHz and 6.30 ± 0.08 /zK R'^ for 
the Finkbeiner Ha template, relative t o 23 GHz. Simi- 
lar res ults were obtained when using the ' Dickinson et al.l 
()2003( ) template with a Ha amplitude of 5.81 ± 
0.05 fiK relative to 23 GHz. The Ha amplitude 



is below what is expected (11.4 /iK R for an electron 
temperature Tp. ~ 8000 K but is consistent with previ- 
ous analyses (jPavies et al.l 120061: [Gold et al.ll2009[ ). The 
monopoles and dipoles are given in Table [2] It is inter- 
esting to see that the monopoles are now all consistent 
with 0, except for Q-band, which appears to still have 
a positive residual monopole of 2.6 ± 0.5 /zK. This may 
be due to the improved foreground modeling while the 
dipoles are of a similar magnitude. This is supported by 
ax^ = 30 1555 ±9 compared to = 301667±9 with just 
one additional degree of freedom. The residual dipoles 
still remain at the « 2 — 3 /iK level and may not be 
totally negligible. 

The CMB temperature power spectra, computed with 
and without the additional Ha template, are shown in 
Fig. [51 The differences are negligible relative to the er- 
ror bars. This tells us, assuming the Ha map is a good 
template for free-free emission, that the single power-law 
model is an adequate model for this particular data set 
and the Kp2 mask. This is justified by the fact that 
the total level of free-free emission at high latitudes is 
actually quite small, relative to the CMB and other fore- 
grounds. The Ha intensity at high latitude is ~ 1 R, 
which corresponds to ~ 10 /xK at 23 GHz and ~ 0.5 fiK 
at 94 GHz. For comparison, at 23 GHz and at high lat- 
itudes, the foregrounds are ~ 100 /iK while at 94 GHz, 
the signal is dominated by CMB at ^ 30 /iK. The free- 
free component is therefore not a dominant component 
at WMAP frequencies, except in the Galactic plane and 
in regions of intense star formation. Indeed, the resulting 
CMB map appears (by eye) to be identical to the single 
template model. However, the difference map does reveal 
that there is a residual Galactic component, as shown in 
Fig. O This residual is ^ 1 /iK for most of the sky, 
increasing nearer the plane to ^ 5 /iK, hence why it is 
not visible in the CMB map itself. The structure of this 
residual component (Fig. [5]) is similar to the Ha map. 
Although it does not have an impact on the CMB power 
spectrum, such errors may have an impact on sensitive 
non-Gaussianity tests. 

The effect on the low-frequency foreground spectral in- 
dex is shown in Fig. [TOl For much of the sky, the differ- 
ences are small (< 0.1), with the largest variations oc- 
curring near well-known regions of strong free-free emis- 
sion such as Orion, Gum nebula and Cygnus. In these 
regions, spectral index decreases by ~ 0.5, as expected 
with a transition from synchrotron to free-free emission. 

4.2.2. Ha free template 

Rather than assume the known spectral dependence 
for a template (the "index" template option), Commander 
allows for an independent fit to each frequency channel 
(the "free" template option) i.e. fitting for a template 
amplitude coefficie nt for each f r equen cy. When this op- 
tion is used for the iFinkbeineil ()2003[ ) Ha template, we 
recover the following amplitudes: 5.67, 2.72, 1.78, 0.68, 
0.27 (±0.08) /iK R-\ at 23, 33, 41, 61 and 94 GHz, 
respectively. These values are slightly lower than those 
obtained with a fixed spectral index (6.3 ± 0.08 /xK R~^ 
at 23 GHz). The amplitudes follow a free-free spectrum 
very closely, with a best-fit spectral index to the am- 
plitudes oi P — —2.12 ± 0.05, consistent with our as- 
sumption of (3ff — —2.15. These coefficients are con- 
sistent with a power-law, with a slight steepening at 
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Fig. 7. — CMB temperature power spectra for the WMAP 5-yr data employing different Galactic masks (percentage values represent the 
sky area masked out). 



1 index template 
3 index templates 




1" 



mmm\m 




10 20 30 40 50 

Multipole f 

Fig. 8. — Top: CMB temperature power spectra for the basic 
WMAP5 model (1 "index" template) and when an additional H« 
template is included (2 "index" templates). Bottom: Difference 
power spectrum between 1 and 2 "index" templates. 

V- and W-ba nds, which is actually expected for free- 
free emission ([Dickinson et alJ l2003f ) . For this model 
we find no evidence for a broad bump, centered at 
^ 40 — 50 GHz in the Ha amp li tudes, t hat has been 
claimed by iDobler &: Finkbeined (|2008bD : 'Doble r et all 
((200i) . This may be due to this particular effect, only 
being visible at lower Galactic latitudes, and hence not 
seen with our larger mask. Indeed, we note the large 
variation in template coefficients when different masks 
are used ( ii4.5|) . 

4.3. Thermal dust 

Thermal (vibrational) dust emission comes from warm 
(Td ~ 18 K) dust emitting black-body radiation, peaking 
at ~ 100 /im (3000 GHz). Thermal dust is a relatively 
weak component in most the WMAP channels. In the 
lowest bands (K,Ka,Q) it is almost negligible, while at 




Fig. 9. — Difference map of the mean CMB maps from the basic 
WMAP5 model (1 template) and when an additional Ha template 
is included (2 templates). Units are fiK. 




Fig. 10. — Difference map of the low frequency spectral indices 
derived from the basic WMAP5 model (1 template) and when an 
additional Ha template is included (2 templates). The full range 
is -0.083 to 1.3. 

W-band it begins to dominate over the other compo- 
nents. One of the assumptions that is commonly used, 
in analyses of WMAP data, is that the thermal dust can 
be subtracted using a single FIR template (e.g. FDS99) 
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- either via a cross-correlation with each band, or with a 
fixed spectral index. This appears to work well, and sur- 
prisingly, the amplitude is remarkably close to the that 
of the template i.e. « 1.0 at 94 GHz (see The 
best-fitting model of FDS99 (model 8) has an effective 
spectral index of /3d « -1-1.7 at WMj4P frequencies. Dust 
is known to emit with a range of indices , with a total 
range of ~ 1 - 3 (jFinkbeiner et al.l Il999t iDupac et all 
120031) but with an average of ~ 2 in the diffuse ISM. We 
now investigate the sensitivity of the results on the as- 
sumption of a spectral index for the thermal dust, and if 
other templates, can provide a better fit. 

4.3.1. The thermal dust spectral index 

To test the sensitivity on the assumed spectral index of 
the dust template, we repeated the WMAP 5-yr analysis 
but varying the assumed dust template spectral index 
with /3d = -t-1.5,1.6, 7, 1.8,1.9,2.0,2.1,2.2, covering the 
likely range for the average value. To reduce the effects 
of free-free emission, which could be sig nificant at the 
higher frequencies, the iFinkbeineil ()2003l ) Ha template 
was included. 

In terms of goodness-of-fit, the total values were 
302240, 301890, 301555, 301234, 300928, 300638, 300364, 
300106 (±440), for increasing /3d. This is teUing us that 
the steeper the dust index, the better the fit. We there- 
fore tried higher values, /3d = 2.5, 3.0, 4.0, 5.0, 6.0... 10.0. 
The minimum of occurs at /3d = -1-6.0. This is be- 
yond the range of plausible physical values (/3d « 1 — 3) 
and must be a consequence of modeling errors. A possi- 
ble explanation is that residuals due to the assumption 
of a single power-law for the low-frequency foregrounds 
(which is accounting for both synchrotron and anoma- 
lous dust emission) are biasing the thermal dust tem- 
plate fit. More precisely, it is known that the anomalous 
emission is tightly correlated wit h thermal dust emission 
based on FIR data ijPavies et al.l [2006) . If the anomalous 
emission is indeed due to spinning dust, the spectrum 
will drop off much more rapidly than synchrotron in the 

40 — 100 GHz range, resulting in a flatter spectral in- 
dex being fitted in the lower frequency channels. This in 
turn will result in too much of the dust-correlated fore- 
ground being subtracted at the higher WMAP frequen- 
cies (V- and W-bands) and thus the dust template will 
prefer a steeper index. The resulting dust spectral index 
iPd = +6.0) can, in principle, give information on the 
residual and hence the spectrum of the modeled compo- 
nent. However, the effective spectral index for the com- 
bined synchrotron and spinning dust components is very 
sensitive on the exact peak of the spinning dust spec- 
trum, since it is close to the lowest frequ ency channels of 
WMA P, in the range ^ 20- 40 G Hz (Dr aine fc LazarianI 
[l998allbl : lAli-Haim oud et al.ll2008ll . Still, given that the 
majority of spinning dust models fall off steeper than 
/3 = —4, it is not inconceivable, that a thermal dust spec- 
trum of /3d — +6.0 is compatible with such a residual. 
The impact on the CMB temperature power spectrum is 
negligible because of the absolute level of thermal dust 
emission at WMAP frequencies is relatively small. 

4.3.2. Thermal dust templates 

So far, we have assumed that the FDS99 template, at 
94 GHz, is the best spatial template of thermal dust at 



WMAP frequenci es. But given the simph city of the fit- 
ted model used bv lFinkbeiner et all ()1999l ) (2 component 
dust model), one can imagine that there may be some er- 
rors in the extrapolation. To see how sensitive the results 
were on this template, we replaced the FDS99 template 
with other templates. First, we tried the E{B — V) red- 
dening maps, which are temperature corrected, and are 
essentially proportional to the dust column density. The 
results were almost identical to the FDS99, with a small 
increase in the value; 295842 for the E{B — V) case 
compared to 295355 in the FDS99 case. This indicates 
that, at least spatially, the E{B — V) map and the FDS99 
maps are very similar. 

We also tried the raw SFD98 100 /im intensity map. 
For this case, = 350863, which is significantly higher 
than in previous cases. Visual inspection of the maps 
revealed that a few pixels, in the Orion region (near 
Barnard's loop), were anomalous. The GMB map was 
strongly negative (—494 /iK) in this region, while the 
low frequency spectral index was at —4.0 i.e. pegged at 
the prior. Also, the power spectrum showed additional 
power at ^ < 10, indicative of foreground residuals. With 
a larger mask, this disappeared. As expected, the FDS99 
template is more representative of the morphology of the 
real sky, than the raw 100 /im intensity map. 

4.3.3. Free dust template 

Up to now, the FDS99 dust template has been used to 
account for the thermal dust contribution at higher fre- 
quencies; we fitted for a template with a fixed spectral 
index e.g. /3 = -1-1.7. We now relax this constraint by 
fitting for the same FDS99 template "freely" at each fre- 
quency. The most important aspect of this model, is that 
the FDS99 template should trace the anomalous emission 
at lower frequencies, as seen in i i4.5l Several combina- 
tions were tried, such as including a Ha template with 
fixed spectral index or not. The results of these were 
similar with no significant variations in the CMB spec- 
trum. For a fixed spectral index Ha template, and the 
Kp2 mask, the Ha amplitude was 8.02 ±0.04 fj,K R"^ at 
23 GHz, similar to what was found in g3](Table[3|). The 
FDS99 amplitude coeflicients were 0.67, 0.08, 0.136, 0.36, 
1.12 (±0.09) nK /iK-i at 23, 33, 41, 61 and 94 GHz, re- 
spectively. The thermal dust contribution appears to be 
robust, with an amplitude at high frequencies remains 
at close to the predicted level from FDS99; this is in 
contrast to the template only fits ( ij4.5| /Table [3]) where 
the values varied considerably with the exact model and 
mask. The dust contribution then decreases substan- 
tially at frequencies immediately below W-band with a 
minimum at 33 GHz. At 23 GHz, the dust-correlated 
contribution then rises quickly, which is indicative of 
spinning dust. However, the value at 33 GHz appears to 
be much lower than expected for typical spinning dust 
models. This is likely to be due to aliasing of power be- 
tween the template and the foreground captured by the 
power-law term in each pixel. Moreover, the FDS99 am- 
plitude at 23 GHz is « 1/10 of the value found with 
template ouly anal yses in and previous analyses 
(|Davies et al.l |2006() . Indeed, simulations showed that 
the dust template coefficients were unreliable, and were 
very sensitive on the details of the spinning dust model. 
We can speculate that only a fraction of the anomalous 
dust signal is truly correlated with the FDS99 template 
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and therefore the rest is accounted for by the low fre- 
quency power-law component. This would also explain 
why a single low frequency power-law appears to be an 
adequate model, at least for WMAP data, for the combi- 
nation of synchrotron and anomalous dust foregrounds. 

4.3.4. Thermal dust as a power-law at each pixel 

Rather than using a fixed spatial template, we relaxed 
this assumption by fitting for a dust amplitude in each 
pixel. This means fitting for two power-laws in each 
pixel, one with a free spectral index (for the low fre- 
quency foregrounds) and one with a fixed index (for the 
thermal dust). WMAP data do not allow fitting for a 
spectral index as well (both in terms of degrees of free- 
dom and frequency coverage), so this has to be fixed for 
the dust power-law e.g. (3d = +1.7. The first attempt 
looked quite reasonable, although the fitted monopoles 
were biased high (4 — 7 /iK). Since the results so far indi- 
cated that residual monopoles/dipoles were small in the 
WMAP 5-yr data, these were fixed at zero. The CMB 
power spectrum was consistent with the previous results. 

The marginal mean map for the dust amplitude is com- 
pared with the FDS99 template in Fig. [TlJ You can see 
that, due to the limited frequency range of WMAP, the 
dust amplitude is relatively noisy at high latitudes. How- 
ever, there is a good overall correspondence between the 
two maps; most of the brighter features are in both maps, 
with similar amplitudes, particularly at lower latitudes 
where the emission is stronger. The FDS99 does a re- 
markable job of predicting the thermal dust amplitude 
at 94 GHz. We did not see a significant improvement 
in the map, except inside the mask (in the Galac- 
tic plane), where the FDS99 prediction is unlikely to be 
accurate due to the complexities of a wider variation in 
dust properties in the Galaxy. Additional high frequency 
channels, as will be provided by the upcoming Planck 
satellite, will provide a much better view of the thermal 
dust component at frequencies > 100 GHz. 

4.4. The spectral index prior 

An important parameter is the Bayesian prior that is 
imposed on the parameters, and in particular, the spec- 
tral index of the power-law component. Commander im- 
plements two types of priors on the spectral index. The 
first is a uniform prior that defines the range of valid 
spectral indices. For the low frequency foregrounds, we 
typically use a range that encompasses any plausible 
value e.g. —4.0 < /? < —1. The second is a more restric- 
tive Gaussian prior, defined by a mean spectral index and 
standard deviation. The advantage of such a prior is that 
in the high signal-to-noise regime, where the data may 
prefer a value several standard deviations away from the 
mean, the solution will move towards the preferred value 
since the Gaussian prior never goes completely to zero 
(until it hits the uniform prior). An example of this can 
be seen in Fig. [3] where near the plane, (3 ~ —2.1, which 
is 3(7 away from the prior mean. On the other hand, in 
the low signal-to-noise regime, the prior becomes more 
important, defining the range of values that the samples 
can take. It is therefore important to choose a prior that 
represents the best available knowledge without being 
too weak or too strong. 
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Fig. 11. — Top: The FDS99 thermal dust prediction template 
at 94 GHz. Bottom: Marginal mean map for the thermal dust 
amplitude (normalized at 94 GHz), when modeled as a power-law 
with fixed spectral index, /3(j = +1.7. 

For the analyses so far, we've used a Gaussian prior^^ 
with mean (3 = —3.0 and standard deviation A/3 — 0.3. 
We now investigate what impact this assumption has on 
the results. 

4.4.1. Mean 

We ran the basic model (as in Q but varying the 
Gaussian prior mean, p. We tried —3.3 < (3 < —2.7 
in steps of 0.1. The total values were 284519, 
284030,283901,284157,284839, 286012, 287771 (±434), 
indicating that /3 = —2.9 gave a marginally better fit. 
The CMB temperature power spectrum remains un- 
changed over this range, showing the robustness of the 
model. 

The FDS99 dust template amplitude varied from 0.978 
(for P = -2.7) to 0.892 (for /3 = -3.3), showing that 
the choice of prior is relatively important for interpret- 
ing the individual foreground components. Similarly, the 
monopoles varied systematically as a function of (3, by up 
to « 10 fiK over the range of priors used here. The cor- 
relation between (3 and the monopole amplitude has neg- 
ligible effect on the CMB power spectrum since it only 
affects the overall offset in the map. 

Fig. [12] shows the foreground spectral index map de- 
rived for P — —2.7, —3.0, —3.3. The bright regions of 
emission do not change by much (A/3 < 0.1), show- 
ing the robustness of the result irrespective of the prior. 
However, at high latitudes, there is a significant differ- 

In practice, the priors need to be multiplied by the appro- 
priate Jeffreys ignorance prior to account for the differen t volumes 
of para meter space for non-linear parameters; see lEriksen et al.l 
Il2008al) for details. 
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ence, indicating that the prior is important in regions 
of weak emission. The most striking difference is in the 
broad ~ 180° diameter "ring" that is very clear for the 
/3 = —2.7 prior, evident for the /? = —3.0 prior, but is not 
visible for the /3 = —3.3 prior (Fig. [T2|) . The ring feature 
is unlikely to be emanating from a particular structure in 
the Galaxy. It is probably a consequence of low signal- 
to-noise regions lying outside the ring, which will tend 
towards the prior mean (e.g. /? = —2.7 or (3 = —3.0) and 
the steeper spectrum regions lying inside, coming from 
the radio spurs (the North Polar Spur) and flatter spec- 
trum emission (e.g. free-free) near the Galactic center. 
When the prior mean is at = — 3.3, this matches closely 
the spectrum of the nearby emission making the feature 
largely disappear (Fig. [T2| . We also note that around the 
Galactic center region, outside the Kp2 mask, the low fre- 
quency spectral index is steeper than average (/? « —3.3), 
which is in disagreement to the hard synchrotron hy- 
pothes is to account for the "haze" (jPobler fc FinkbeineH 
12008a') . A similar analysis was repeated for the case with 
2 templates, and similar results were found. 

4.4.2. Standard deviation 

We now vary the standard deviation, keeping the mean 
fixed at /3 = -3.0. We tried A/3 = 0.1, 0.2. ..0.5 in steps 
of 0.1. The values were 313191, 289437, 284157, 
282694, 282280 (±434). As expected, the goodness-of-fit 
increases steadily with increasing A/?, but with a signif- 
icant jump for A(3 < 0.2. This suggests that the prior is 
dominating for A/3 < 0.2 and hence there are real fluctu- 
ations of A/3 > 0.2 on the sky. The maps of the spectral 
index mean for Af3 — 0.1,0.3,0.5 are shown in Fig. [T3l 
It is clear that A/3 = 0.1 strongly dampens the variation 
of spectral index at high latitudes. Our original choice 
of A/3 = 0.3 is justified since there is only a very small 
change from A/3 = 0.3 to A/3 = 0.5. 

4.4.3. Running of the spectral index 

The assumption of a power-law, for the low frequency 
foregrounds, is unlikely to be accurate over a wide range 
of frequencies. Synchrotron emission is known to steepen 
with frequency due to spectral aging of the cosmic 
ray electrons, although multiple synchrotron components 
(with different spectral indices) may result in a flatten- 
ing of the effective index. Anomalous emission could also 
modify the effective spectral index, probably resulting in 
a steepening at higher frequencies. A first order correc- 
tion to the power-law is to fit for an additional parameter 
in the exponent, usually referred to as the "running of the 
index", or curvature parameter, C (see Tabic [ij. Note 
that positive values of C represent a flattening of the 
spectral index with increasing frequency, while negative 
values represent a steepening with frequency. 

The WMAP data alone do not allow for a reliable 
estimate of C for each pixel. Instead, we fit a model 
with a fixed value of C for the whole sky, relative to 
23 GHz. We found that there was a slight preference for 
a flattening of the spectral index; the values where 
295362, 293725, 292757, 292386, 292776, 294427 (±450), 
for C = 0.0,0.1,0.2,0.3,0.4,0.5, respectively. The best- 
fitting value therefore appears to be for C — 0.3 i.e. 
a flattening of the spectral index at higher frequencies. 
This is also the preferred value found by iKogut et al.l 
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Fig. 12. — Low frequency foreground spectral index maps 
for different Gaussian prior means. From top to bottom, /3 = 
—2.7, —3.0, —3.3 with a standard deviation, A/3 = 0.3. 



(|2007f ) for the polarized synchrotron, although they 
found C « —0.4 for the unpolarized emission. The dis- 
crepancy is likely to be due to modeling errors and the 
additional components in temperature, that do not emit 
strongly in polarization (e.g. free- free, spinning dust). 
The flattening of the index is not what one would initially 
expect in temperature, since a strong spinning dust com- 
ponent would likely result in a negative curvature. There 
are several explanations that could account for this: i) 
multiple synchrotron components with a range of spec- 
tral indices would result in a flattening of the index at 
higher frequencies, ii) free-free emission unaccounted for 
by the Ha template, iii) thermal dust emission not ac- 
counted for by the FDS99 templat e . It is interesting 
to note that Ide Oliveira-Costa et al.l ()2008f ) observed an 
overall flattening (C > 0) when fitting data from 10 MHz 
to 100 GHz, adopting a reference frequency of 5 GHz. 
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Fig. 13. — Low frequency foreground spectral index maps for 
different Gaussian prior standard deviations. From top to bottom, 
A/3 = 0.1, 0.3, 0.5 with a mean, 13 = -3.0. 

4.5. Template only fits 

The fitting of spatial templates to CMB data has 
been one of the most effective wa ys of reducing a.nd un - 
derstanding foreground emissions (|Hinshaw et al.ll2007D . 
The template amplitude coefficients can also give addi- 
tional insight into the scaling of the foregrounds. We 
investigate the performance and results of fitting tem- 
plates only to the WMAP 5-yr data. 

4.5.1. Haslam 408 MHz, Ha and FDS99 

We be gin by using " s tanda rd" templates for 3 com- 
ponents: iHasiamleLan ()1982f ) 408 MHz as a tracer of 
synchrotron, Ha for free-free emission and FDS99 for 
thermal dust (the FDS99 template also accounts for 
the low frequency anomalous dust-correlated emission). 
Each template is fitted to each frequency channel in- 
dependently, using pixels outside the Galactic mask. 
Monopoles and dipoles were not fitted for. 



The template amplitude coefficients are given in Ta- 
ble [3] for different Galactic masks. It can be seen that 
there are significant variations in the coefficients, partic- 
ularly at the higher frequencies. It is clear that strong 
emission just outside the mask can have a strong ef- 
fect on these single-fit numbers. Overall, the num- 
bers are similar to those found by pre v ious authors 
(iBa ndav et all l2003t iBennett et all 120031 : iDavies et all 
2006). The FDS99-correlated foreground is a minimum 
at V-band but increases quickly at lower frequencies. 
This result is one of the key arguments for an anomalous 
component while the spinning dust origin naturally ex- 
plains the strong correlation with traditional dust emis- 
sion, as traced by the FDS99 template. The spectral 
index of this foreground, in the WMAP frequency range, 
is /? w —4 which is consistent with a spinning dust type 
of spectrum. The individual coefficients correspond to 
a synchrotron spectral index of /3s ~ —3.0, consistent 
with a slight steepening relative to GHz frequencies. The 
408 MHz coefficients on the other hand fall off somewhat 
slower than what would be expected from synchrotron 
emission, with /3s « —2.5, probably due to aliasing of 
power between the templates and the low signal of syn- 
chrotron emission at the higher WMAP frequencies. The 
Ha coefficients also fall off slower than expected, with 
/3 ~ —1.4, rather than the expected value for free-free 
emission, f3ff = —2.15. Again, this is likely to be due 
correlations between the templates. Still, the Ha coeffi- 
cients are significantly lower than what is expected from 
theory, by almost a factor of 2. For example, at 22.8 GHz 
we obtain coefficients in the range 6.67 — 8.52 /iK 
at 22.8 GHz, while for a typical electron temperature 
of Te = 8000 K the theoretical value is 11.4 fiK R''^ 
(jPickinson et al.ll2003f ). The origin of this is still not un- 
derstood, but is likely to be related to the correlations 
between the various components. 

The CMB power spectrum from this template analysis 
and for the Kp2 mask is shown in Fig. 1141 The power 
spectrum is clearly strong contaminated by foreground 
emission with a significant positive bias, particularly at 
£ < 20. This is largely because the synchrotron emission 
at 23 GHz is not perfectly traced by the 408 MHz map. 
Spatial variations in the synchrotron spectral index will 
effect the morphology of the map resulting in an imper- 
fect template. Moreover, the 408 MHz map is known to 
have baseline offsets resulting in stripes in the final maps. 
Finally, there are residual extragalactic sources that have 
not been accounted for, although at the working resolu- 
tion of 3° the impact of these residuals is expected to be 
small. 

4.5.2. K-Ka template 

As shown above, using the "standard" templates does 
not fully remove all the foreground emissions on large 
angular scales. A better synchrotron template, pro- 
posed by the WMAP team, is the K— Ka difference map 
(|Hinshaw et al.|[2007f ). This has the advantage of per- 
fectly canceling out the CMB and thus containing only 
foreground emission at an effective frequency within the 
WMAP frequency range. Furthermore, the absolute cal- 
ibration will be perfectly matched to the other WMAP 
frequencies. The disadvantages are i) the K-band and 
Ka-band maps cannot be used for cosmological analyses 
and ii) the template amplitude coefficients are more diffi- 
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TABLE 3 

Template amplitude coefficients for a 3 template fit of Haslam et al. 408 MHz, Finkbeiner Ha and FDS99 94 GHz. Units 

ARE K K^^, fiK R"^ AND /.tK {fj,K)~^, respectively. 



Template Band 


Mask 


408 MHz K 
Ka 

Q 
V 

w 


4.34 ±0.08 4.24 ±0.11 4.43 ±0.13 4.57 ±0.08 
1.67 ±0.08 1.80 ±0.11 1.88 ±0.13 1.81 ±0.08 
1.02 ±0.08 1.20 ±0.11 1.26 ±0.13 1.15 ±0.08 
0.42 ±0.08 0.64 ±0.11 0.69 ±0.13 0.55 ± 0.08 
u.zo in Kj.vio U.10 in Kj.A-± u.oz zn u.±o u.ot u.uo 


Ha K 

Ka 

Q 
V 

w 


6.67 ±0.09 6.76 ±0.16 7.38 ±0.14 8.52 ± 0.09 
3.87 ±0.09 4.62 ±0.16 4.49 ±0.14 5.20 ± 0.09 
2.82 ±0.09 3.67 ±0.16 3.37 ±0.14 4.20 ± 0.09 
1.62 ±0.09 2.50 ±0.16 2.11 ±0.14 2.74 ± 0.09 
1.19 ±0.09 2.06 ±0.16 1.64 ±0.14 2.26 ± 0.09 


FDS99 94GHz K 
Ka 

Q 
V 
W 


6.03 ±0.11 5.64 ±0.11 5.58 ±0.16 5.62 ± 0.09 
1.79 ±0.11 1.15 ±0.11 1.23 ±0.16 1.33 ± 0.09 
0.81 ±0.11 0.11 ±0.11 0.22 ±0.16 0.30 ± 0.09 
0.36 ±0.11 -0.36 ±0.11 -0.25 ±0.16 -0.18 ±0.09 
0.89 ±0.11 0.13 ±0.11 0.26 ±0.16 0.34 ± 0.09 




Multipole £ 

Fig. 14. — CMB temperature power spectra from Commander in- 
cluding the basic model, 3 traditional templates (408 MHz, 
Ha, FDS99 94 GHz) and 3 WMAP templates (K-Ka, Ha, 
FDS99 94 GHz). All analyses are using the KQ85 temperature 
mask. 

cult to quantify, since they are effectively a contribution 
of at least two (or more) Galactic components. 

We analyzed the WMAP5 data using the K— Ka dif- 
ference template, fixing the dust spectrum to /3d — +1.7 
and the free-free spectrum to (3g — —2.15. The cor- 
responding CMB power spectrum, for the KQ85 mask, 
is compared with the traditional templates in Fig. [Ml 
The power spectrum is very close to the basic WMAP5 
analysis (fj3]) indicating the improvement in using the 
K-Ka map over the Haslam et al. (1982) 408 MHz. 
The improvement is partly due to the more accurate low 
frequency template, but also the fact that the K- and 
Ka-bands can no longer be used because these data are 
the bases of the template. The template amplitude co- 
efficients are given in Ta ble [H The numbers are close 
to the values from lHinshaw et al.l ()2007f ) but not 

exactly the same. The differences are likely to be due to 
the lower resolution of our analysis (3° compared to 1°) 
and the modification of the temperature mask due to the 
additional smoothing used here. We also tried fitting of 
monopoles and dipoles with similar results. 

4.6. CMB pixel fitting 

The results from the previous sections included full-sky 
foreground amplitude and spectral index maps where the 
low frequency foreground is fitted inside the mask. One 
issue when interpreting the foreground results is that 



TABLE 4 

Template amplitude coefficients for an analysis with Q-, 

V- AND W-BANDS ADOPTING THE K-KA DIFFERENCE MAP AS A 
LOW FREQUENCY TEMPLATE. UNITS ARE K , flK AND 

^J.K (/iK)-i, FOR THE K-Ka, Ho AND FDS99 94 GHz 

TEMPLATES, RESPECTIVELY. THE KQ85 MASK WAS ADOPTED. 



Band 


K-Ka 


Ha 


FDS99 94 GHz 


Q 


327 ± 14 


1.21 ±0.11 


1.21 ±0.11 


V 


115 ± 14 


0.60 ±0.11 


0.60 ±0.11 


w 


50 ±14 


0.35 ±0.11 


0.35 ±0.11 



only the large scale CMB fluctuations are constrained 
within the mask, resulting in CMB residuals inside the 
masked area that can effect the foreground reconstruc- 
tion. This is the reason for the discontinuity between 
the masked and unmasked region, as seen in Fig. [3] To 
obtain a more accurate foreground map. Commander can 
also fit for the CMB in the pixel domain. In this model 
the CMB is simply fitted by an amplitude in each pixel, 
since the CMB spectrum is known (Table [it; this is sim- 
ilar to the MCMC algorithm described by lErik s en et al.l 
(2006). Each pixel is treated as independent, therefore 
spatial correlations are not taken into account, which can 
result in a loss of information and hence noisy solution 
at small angular scales. 

Fig. [T5| shows the margin mean maps for the CMB, 
foreground amplitude and index, and their associated 
r.m.s. uncertainties. For most of the sky, the results are 
very similar to the results already presented. However, 
there are significant differences inside the mask due to 
the different treatment of the CMB in these pixels. The 
CMB map reveals where the model/frequency range is 
not adequate. First, the strong free-free regions are seen 
as bright pixels in the CMB map (e.g. Cygnus at I ~ 90°) 
due to the flatter free-free spectral index. In the Galac- 
tic center region, there is a slight lack of power. Exam- 
ination of the spectrum at these pixels reveals nothing 
special, except for a flattening at W-band due to ther- 
mal dust. The likely reason for this is the mis-scaling of 
the FDS99 94 GHz dust template at low latitudes. In 
the Galactic plane, there is a complicated mix of dust 
properties (composition, temperatures, emissivities etc.) 
which is unlikely to be well-represented by the simple 
2-component fit made bv iFinkbeiner "eTaLl dHH). 
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Rather than using a fixed spatial template, it is also 
possible to fit for the thermal dust component as an ad- 
ditional foreground to be fitted in each pixel. Clearly 
there are not enough channels to fit for a more compli- 
cated model. What can be done, is to fit for the ther- 
mal dust amplitude with a fixed spectral index (as in 
t j4.3.4p . This reduced the negative values observed near 
the Galactic Center. We found that /9d = 2.0 gave a 
better fit with =1804915 co mpare d to = 1885888 
for /3d = +1-7. As we saw in ! j4.3.11 but in a different 
context, the preferred value is actually steeper than this 
and reflects the unmodeled flatter spectral index (e.g. 
free-free) emission that is biasing the low frequency spec- 
tral index. The monopoles and dipoles for this model 
are in good agreement with the previous models, al- 
though they are statistically much more accurate (ig- 
noring modeling errors). From low to high frequency 
they are: 0.177 ± 0.003, -1.844 ± 0.011, 3.188 ± 0.007, 
-1.930 ± 0.011 and 0.408 ± 0.003 /^K. 

5. COMPARISON OF CMB AND FOREGROUNDS WITH 
OTHER WORKS 

5.1. CMB maps and low-l power spectra 

In this paper, we have produced several versions of the 
CMB map with different assumptions on the foreground 
signal model. We now make a more detailed comparison 
between these maps and those available in the literature. 
We take our CMB map from the 2 template plus power- 
law model f ^4.2.1l) and comp a re it with i) the WMAP 
5-yr ILC map Hinshaw ct al. (2009*), ii) the Harmonic 
ILC map of [Kim ct al. (2008) and ii i) the WMAP 5-y r 
MCMC-derived best-fit CMB map of lGold et all ((20091 ). 
The maps were smoothed to a common 3° resolution and 
downgraded to A'sidc — 64. The MCMC map is already 
at iVside = 64 thus we simply smooth to 3° resolution. 

A by eye visual inspection of the CMB temperature 
maps reveals no obvious differences, except in the Galac- 
tic plane, which is clearly strongly contaminated by fore- 
grounds. To make a more accurate comparison, we mask 
the Galactic plane and make difference maps as shown 
in Fig. 1161 

There is good overall agreement with the WMAP ILC 
(WILC) map with very little large scale power. There are 
a few regions close to the Galactic plane mask and LMC 
that are clearly different, with more power in the ILC 
map (blue regions on Fig. [TO)) . This to be expected due 
to the extra degrees of freedom afforded by the power- 
law at each pixel, compared to the weights of the ILC 
method which are computed over large regions. The most 
apparent difference is the significant amount of small- 
scale power, which is a consequence of the bias correction 
term employed in the WILC map, resulting in significant 
"noise" at scales of ~ 1°. 

However, there appears to be significant differences, 
even on large angular scales, between t he Command e r map 
and the Harmonic ILC map. Indeed, iKim et al] ()2008D 
noted that most of the low-^ powers of the HILC map are 
higher than those of the WILC map. This suggests that 
the HILC map contains some residual foreground power, 
as seen in Fig. [1^1 

The WMAP MCMC-derived map appears to show a 
significant positive bias. Outside the mask, the MCMC 
map has a mean of —8.6 /iK while the difference map has 
a mean of -f8.5 /xK. Furthermore, there are large scale 



features (particularly in regions of known strong Galactic 
emission) that can be seen in the difference map, that 
are likely to be due to foreground modeling errors and 
residual mo nopoles/dipoles. 

Recently, iPietrobon et al.l ()2008f ) used needlets to de- 
tect non-isotropic features in the CMB sky i.e. non- 
Gaussianity. By masking out the brightest (3(t) detec- 
tions of hot and cold spots, they claimed that a sig- 
nificant reduction in the CMB power occurred between 
£ = 8 — 30 at the \(j level. We repeated our analysis 
using the additional areas masked out by the hot/cold 
spots. The CMB temperature power spectra are com- 
pared in Fig. [T71 The spectra are in good agreement 
with the standard KQ85 mask, except at ^ < 20 where 
we confirm a slight decrease (« 100 /iK^) in power. 

5.2. North vs South ecliptic asymmetry 

A number of analyses of WMAP data have revealed de- 
partures from statistical isotropy on large angular scales. 
One of the most apparent, is the CMB power asymme- 
try between the northern and southern ecliptic hemi- 
spheres. Besides for possible low-l alignments, there ap- 
pears to be more power in the so uthern hern i spher e, than 
in the north (see for exam ple, Ha nsen et al.l ()2004l . [20081: 
lEriksen et all (l200l . [2007i) : iBernuil (|2008')). 

We ran Commander with the same basic model as used 
in Sj3l using a number of different masks for the eclip- 
tic areas, based on a cut at ecliptic latitude 30° e.g. 
the North ecliptic Kp2 covers ecliptic latitudes > 30° 
in addition to the Galactic Kp2 mask, while the eclip- 
tic plane is for latitudes within 30° of the ecliptic plane. 
The left plot in Fig. [T51 shows the Commander CMB map, 
with Kp2 mask in grey, and ecliptic region boundaries 
shown as black lines. We also tried the same using the 
WMAP KQ85 temperature mask and the 10° smoothed 
Kp2 mask (see i )4.ip . 

Fig. [TH] shows the power spectra for a selection of 
masks, based on the Galaxy and ecliptic latitudes. For 
clarity, we plot binned power spectra with A£ = 5. The 
North-South asymmetry is visible, mostly at £ < 15. The 
three lowest £-bins are « 2a discrepant between the north 
and south ecliptic hemispheres. The CMB in the eclip- 
tic plane is in between, but closer to the southern val- 
ues. Interestingly, the power spectra for the southern 
ecliptic hemisphere, when using a larger Kp2 smoothed 
Galactic mask contains less power at ^ < 10. This can 
be understood by simply looking at the CMB map and 
observing where the large-scale features lie relative to 
the ecliptic (left hand of Fig. [T5|) . The two bright large- 
scale features are both near the edge of the standard Kp2 
mask; a hot spot near south of the Gum nebula region 
((/, 6) w (260°, -20°)) and a cold spot to the south-west 
of the Galactic center, at ecliptic latitude « —30°. Both 
of these will be largely masked out in the smoothed Kp2 
mask, thus reducing the power at low-i. However, both 
of these features appear to remain robust against fore- 
ground subtraction and therefore remain a mystery. 

5.3. Synchrotron / low frequency foreground spectral 
index maps 

One of the products of running CommEinder, besides 
the CMB posteriors, are the posteriors for the fore- 
ground parameters. In 52] we presented the CMB and 
foreground results when applied to WMAP 5-yr data. 
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Fig. 15. — Reconstructed mean and r.m.s. maps for the CMB, low frequency foreground amplitude and spectral index. Here, no mask 
was used and the CMB is fitted pixel-by-pixel with an amplitude. 



An important product is the spectral index of the low 
frequency foregrounds, as fitted over the WMAP fre- 
quency range, for each pixel. This accounts for all the 
diffuse foreground components at low frequencies, and 
could include contributions from synchrotron, free-free 
and anomalous emission. There now exists several ver- 
sions of synchrotron/low-frequency spectral index maps 
available in electronic form. 

Fig. [19] shows spectral index maps from several au- 
thors, smoothed to a common resolution of ~ 3°. Table[5] 
lists the maps, the frequency range, reference and the 
mean/r.m.s. spectral index at high latitude {\b\ > 20°). 
The top row of maps in Fi g.llOlwere derived u sing low fre- 
quency data at 408 M Hz (|Haslam et al.lll982D. 1420 MH z 
(|R,eich fc R,eichl[T98fi[ ) and 2326 MHz (|,Tonas et all flQQSh. 
The bottom row includes WMAP data, either on their 
own, or combined with low frequency data (Table [5]) . 

There is considerable variation between these maps, 
which can be partly attributed to the different frequency 
ranges used, and also different assumptions that were 
made in their production. The mean spectral index and 
r.m.s. variation at high latitude > 20°) are given in 
Table [S] The low frequency spectral index maps are typ- 
ically flatter, with (3 w -2.8, while at WMAP frequen- 
cies, P ~ —3. 0. This is in br oad agreement with pre - 
vious studies (iReich fc R eich 19881 iDavies et al.lll996h . 
Ide Qliveira-Costa et al.l ((20081) analyzed radio maps from 
10 MHz to 100 GHz and found a -2.6 < f3 < -2.3 at 
150 MHz, but -3.3 < l3 < -2.1 at 5 GHz. This is con- 
sistent with what is expected for synchrotron radiation, 
where spectral aging of cosmic ray electrons reduces the 
power at higher frequencies. 

The low frequenc y spectral index maps of Finkbeiner^^ 
(priv. comm.) and iPlatania et all (|2003f ) are very simi- 
lar. They used similar datasets and Fourier filtering tech- 
niques for reduci ng the striping artifac ts that exist in the 
radio maps. The lGiardino et al.l (|2002 ^ map is somewhat 
different, both in absolute terms but also in morphology, 
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due to different processing techniques, particularly for 
dealin g with str iations in the radio maps. 

The lGold et al. ( 2009 ) spectral index map has a similar 
morphology to the Commcinder map, with steeper indices 
towards the Galactic center region, and flatter indices 
outside. The Commander map is typically flatter in the 
Galactic plane region, which is probably a result of using 
spatial templates for fr ee-free/thermal-du st, compared to 
the MCMC analysis of iGold et all ([2009h . The Giardino 
(priv. comm.) 408 MHz-23 GHz spectral index map is 
somewhat different in nature. Much of this is attributed 
to the larger frequency range employed a nd the reliance 
on th e WMAP synchrotron MEM model (|Hinshaw et all 
I2OOI as a reliable map of the synchrotron radiation 

We note once again that we do not see a signifi- 
ca nt evidence f or a " WMAP Haze", as was first seen 
by iFinkbeiiie^ ()2004al) . and later suggested to be flat 
spectrum synchrotron from dark matter anni hilation 
(|Finkbeinerj I2004bl: iDobler fc Finkbeineil l2008aD . First 
of all, we do not see a clear flattening of the spectral in- 
dex near the Galactic center region (e.g. Fig.[3l); in fact, 
the integrated foreground is typically steeper in this re- 
gion outside the mask, although the "haze" is only a 
sub-component of the total emission. Secondly, there is 
no clear sign of an extra unmodeled component in the 
goodness-of-fit map of Fig. although we do note 
that the does increase slightly towards the Galactic 
plane (as expected from the increasing complexity) and 
also in the ecliptic plane (due to additional unmodeled 
instrumental noise) which happens to pass through the 
Galactic center region. 

6. CONCLUSIONS 

The Gibbs sampling algorithm provides a self- 
contained method to perform component separation and 
CMB power spectrum estimation for multi-frequency 
data. As implemented in the Commander code, we can 
fit for a variety of signal models, including spatial tem- 
plates and/or parametric functions on a per pixel ba- 

ftp://ftp.rssd.esa.int/pub/synchrotron/README.html 
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Fig. 16. — CMB temperature difference maps, relative to the 
Commander 2 index template plus power-law model. All maps have 
been smoothed to 3° resolution and are displayed with a linear 
scale ranging from -15 /iK to +15 fjK. Top: WMAP 5-yr ILC 
map; Middle: Harmonic ILC map; Bottom: WMAP 5-yr MCMC 
map. 



TABLE 5 

Characteristics of the synchrotron / low frequency 
spectral index maps in fig. [191 the last column gives the 
mean spectral index (/3), and the r.m.s. variation (a/3) at 
high latitude > 20°). 



Frequency 
range 



Reference 



/3± A/3 
(|b| > 20°) 



0.408 -2.3 GHz Finkbeiner (priv. comm.) -2.69 ±0.12 
0.408 - 2.3 GHz Giardino et al. (2002) -2.90 ± 0.06 
0.408 - 2.3 GHz Platania et al. (2003) -2.68 ± 0.11 

23 - 94 GHz Gold et al. (2008) -3.18 ± 0.17 

0.408 - 23 GHz Giardino (priv. comm.) -2.89 ± 0.13 

23 - 94 GHz This paper -2.97 ±0.21 

sis, and propagate errors seamlessly the CMB likelihood. 
The method and code has now been tested and validated 
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Fig. 17. — Top: Comparison of the CMB temperature power 
spectrum when using the KQ85 mas k, and KQ85 m ask plus the 
additional hot/cold spot detections of lPietrobon et all |[2008i ). Bot- 
tom: Difference power spectrum. 

on a variety of simulations and real data. 

In this paper, we applied the Commander code to the 5- 
yr WMAP data using a signal model consisting of CMB, 
noise, a power-law with amplitude and spectral index in 
each iVside = 64 pixel, a spatial template for thermal 
dust, and monopoles/dipoles at each frequency. For the 
Kp2 mask, the recovered power spectrum, up to £ = 50, 
is in very good agreement with the published WMAP 
team spectrum and the cosmological parameters are neg- 
ligibly affected. Using larger masks did not improve the 
power spectrum, in terms of reducing foreground resid- 
uals. However, smaller masks (< 15% masked out) did 
result in a bias in the power spectrum, indicating signif- 
icant residuals. A full-sky analysis results in ^ 500 fiK? 
of additional power for ^ < 50. 

Residual monopoles and dipoles, that were originally 
found to be a problem with WMAP 3-yr data, are found 
to be small in the 5-yr data, within a few /iK of zero. 
Although very small, not accounting for them, can af- 
fect the details of the recovered foreground parame- 
ters. The thermal dust template amplitude at 94 GHz 
is remark ably clo se (within « 5%) to the prediction by 
iFinkbe inc r et al.l ([l999), assuming a dust spectral index 
/3dust = +1.7. A power-law, with amplitude and spec- 
tral index in each pixel, accounts for all the low fre- 
quency diffuse foregrounds, and is found to be an ad- 
equate model outside the Kp2 mask at WMAP frequen- 
cies. The best-fitting thermal dust spectral index, based 
on the goodness-of-fit statistic, appears to be much 
steeper than would be expected, and is likely to be an 
artifact of modeling errors due to dust-correlated emis- 
sion with a different spectrum. Again, the effect on the 
CMB power spectrum is negligible, due to the low abso- 
lute levels of thermal dust at WMAP frequencies. 

We tested the assumptions that were made about 
the signal model including the spectral index prior 
(mean and width), thermal dust spectral index, 
monopoles/dipoles. In all cases, the impact on the CMB 
power spectrum was very small to negligible. However, 
the assumptions did change the details of the recovered 
foreground parameters, particularly at very high lati- 
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Fig. 18. — Left: Commander CMB map, with the Kp2 mask region in grey. The echptic regions (defined by ochptic latitudes above or 
below 1 30° I) are marked. Right: Commander power spectra for the standard Kp2 mask, and various ecliptic regions, combined with Galactic 
masks. Each data point is the power spectrum in t-hins with h.i = 5. Symbols are shifted in t for clarity. 
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Fig. 19. — Comparison of synchrotron / low frequency spectral index maps. Top row: low frequency spectral index maps, from left to 
right: Finkbeiner 408 MHz to 2.3 GHz, Giardino et al. (2002) 408 MHz to 2.3 GHz, Platania et al. (2003) 408 MHz to 2.3 GHz. Bottom 
row: spectral index maps including WMAP data, from left to right: Gold et al. (2009) 23 - 94 GHz, Giardino 408 MHz to 23 GHz, 
Commander 23-94 GHz. 



tudes, where the signal-to-noise ratio in each pixel is rela- 
tively low, and therefore the priors are having an impact. 

To account for free-free emission, we fitted for a Ha 
template and found amplitudes that were lower than the- 
ory (assuming Tg « 8000 K), but is similar to those found 
by other authors. The effect on the CMB map is small 
(< 1 /.tK) for most of the sky but rises to ^ 5 — 10 //K 
at lower latitudes near the edge of the Kp2 mask. Simi- 
larly, the effect on the map of low frequency spectral in- 
dex is also small (A/3 < 0.1) except near well-known star 
forming regions (e.g. Orion, Gum nebula) where free- 
free emission dominates at WMAP frequencies. This is 
as expected, since a single power-law cannot account for 
two components (e.g. synchrotron and free-free) that are 
approximated by power-laws with different indices. 

When fitting for templates only (i.e. no spatial varia- 
tion in spectral index), we find that using a combination 
of 408 MHz, Ha and FDS99 templates, results in resid- 
ual foreground power at low £. This indicates that this 
model is not sufficient to account for all the foreground 
complexities, particularly at the lower WMAP frequency 



channels. We found that the template amplitudes were 
sensitive to the details of the Galactic mask. The FDS99 
template amplitudes rise at the lower WMAP channels, 
indicative of spinning dust, and similar to previous anal- 
yses. Replacing the 408 MHz template with the K— Ka- 
band difference map, and fitting to Q-, V-, W-bands, 
results in a cleaner separation and unbiased power spec- 
trum. 

The Commander CMB map was compared with maps 
available in the literature. Differences of ~ 10 — 20 ^xK 
were observed, and can be attributed to the differences 
in techniques and assumptions of component separa- 
tion. The asymmetry of power between the north and 
south ecliptic hemispheres is confirmed and remains ro- 
bust against foreground modeling. The Commander low 
frequency spectral index has a similar morpholog y and 
range to the MCMC analysis of iGold et all ()2009[ ). with 
13 = -2.97 ± 0.21 at high latitude. The spectral index 
is typically steeper by A/? « 0.2 than is observed at low 
1 GHz) frequencies. 

Commcinder will be a valuable \ow-(. tool for interpret- 
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ing future CMB data, such as the upcoming Planck mis- 
sion. The increased sensitivity, and more importantly, 
larger frequency range will allow a more complete fore- 
ground model to be fitt ed for, and should avoid the need 
for spatial templates (iLeach et alJ 12008). In a forth- 
coming paper, we will investigate the effectiveness of 
WMAP+Planck data when fitting for CMB and multi- 
ple foreground components. The increased sensitivity in 
polarization will allow the same procedure to be applied 
to the Q and U polarization maps. 
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